Coulomb drag between massless and massive fermions. 
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We theoretically investigate the frictional drag induced by the Coulomb interaction between spa- 
tially separated massless and massive fermions in the Boltzmann regime and at low temperatures. 
As a model system, we use a double-layer structure composed of a two-dimensional electron gas 
(2DEG) and a n-doped graphene layer. We analyze this system numerically and also present ana- 
lytical formulae for the drag resistivity in the limit of large and small interlayer separation. Both, the 
temperature and density dependence are investigated and compared to 2DEG-2DEG and graphene- 
graphene double-layer structures. Whereas the density dependence of the transresistivity for small 
interlayer separation differs already in the leading order for each of those three structures, we find 
the leading order contribution of the density dependence in the large interlayer separation limit to 
exhibit the same density dependence in each case. In order to distinguish between the different 
systems in the large interlayer separation limit, we also investigate the subleading contribution to 
the transresistivity. Furthermore, we study the Coulomb drag in a double-layer structure consisting 
of n-doped bilayer and monolayer graphene, which we find to possess the same qualitative behavior 
as the 2DEG-graphene system. 
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I. INTRODUCTION 

The transport properties of double-layer systems, in 
which carriers arc confined to nearby parallel planes, 
have received considerable attention since the earlier pro- 
posal by Pogrebinski! 1 of employing a bilayer system for 
measuring the frictional drag. Drag measurements are 
performed by driving a current j a through one of the 
layers (the active layer) and measuring the electric field 
E p induced in the other layer (the passive layer) due to 
interlayer momentum transfer. The drag transresistiv- 
ity (also called drag coefficient, or simply, drag) is de- 
fined as pa = Ep/j a . The measurement of the frictional 
drag can provide valuable information about the density 
and temperature dependence of the carrier-carrier inter- 
action in 2D systems. In particular, the frictional drag 
due to the interlayer carrier-carrier Coulomb interaction 
in double-layer semiconductor systems has been investi- 
gated in great detail. 2-10 

With the recent progress in the physics of graphene, 
much attention has been devoted to the investigation 
and understanding of the frictional drag in spatially sep- 
arated double-layer graphene systems 11-19 as well as in 
structures comprising two bilayer graphene (BLG) sheets 
isolated from each other by a spacer. 16 Moreover, in the 
limit of low temperatures and large interlayer distances, 
a generic formula has been derived for the leading order 
of the asymptotic behavior in the limit of large interlayer 
separation for systems where each layer I is described by 
an energy dispersion of the form e l w cx k 2 ~^ 1 (£, is a layer 
specific constant) and a momentum-dependent relaxation 
time Ti (k). 19 

There is a fundamental difference between the carriers 
in graphene and those in a two-dimensional electron gas 
(2DEG) or in BLG. While in graphene the carriers can be 



interpreted as massless fermions with a linear dispersion, 
in the 2DEG and BLG the carriers exhibit a parabolic 
dispersion and have a finite effective mass. Thus, most 
of the previously reported investigations of the frictional 
drag have been limited to the case of interaction between 
massive fermions (in the case of 2DEG-2DEG or BLG- 
BLG double-layer systems) or between massless fermions 
(in the case of graphene-graphene double-layer struc- 
tures). In what follows, we will refer to the former and 
later cases as massive-massive and massless-massless sys- 
tems, respectively. By assembling a double-layer struc- 
ture consisting of a graphene layer and a 2DEG layer, it 
might also be possible to create a setup where the car- 
rier densities are significantly different in both layers, a 
case difficult to achieve if both layers consist of the same 
material. 

In the present paper, we investigate massless-massive 
systems in which the frictional drag is induced by 
the Coulomb interaction between massless and massive 
fermions, a case that till now has remained largely un- 
explored. Here, we restrict ourselves to the discussion 
of low temperatures and the case where both layers are 
in the ballistic/Boltzmann regime. As a prototype sys- 
tem, we consider first a double-layer structure consist- 
ing of a 2DEG formed in a GaAs quantum well and a 
closely located n-doped graphene layer. We then com- 
pute the transresistivity for such a system and investi- 
gate its dependence on temperature and carrier concen- 
trations. We also provide analytical formulae describ- 
ing the asymptotic behavior of the transresistivity in the 
large and small interlayer separation limits and compare 
our results with those corresponding to 2DEG-2DEG, 
graphene-graphene, and BLG-BLG double-layer struc- 
tures. We show that in the small interlayer separation 
limit, already at the leading order the transresistivity 
scales with the carrier densities differently for all the 
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FIG. 1: (Color online) Schematic illustration of the geometry 
considered. The 2DEG is located within a quantum well of 
width w and its localization along the z direction is described 
by X2d(z), whereas the location of the graphene sheet is given 
by Xg( z )- The relative dielectric constants of the structure are 
given by Ki, «2, K3, and k 2D . 



three massive-massive, massless-masslcss, and masslcss- 
massive systems. However, in the large inter layer sepa- 
ration limit the three kinds of systems exhibit the same 
asymptotic behavior in the leading order and differences 
appear only when the subleading correction is taken into 
account. As an alternative to the 2DEG-graphene struc- 
ture, we also investigate the drag transresistivity in a 
double-layer structure consisting of n-doped bilayer and 
monolayer graphene isolated from each other by a spacer. 
Such a masslcss-massive system exhibits the same quali- 
tative behavior as the 2DEG-graphene structure. 

The manuscript is organized as follows: In Sec. II, fol- 
lowing the introduction of the model and the theoreti- 
cal framework, the Coulomb drag in 2DEG / (monolayer) 
graphene systems is discussed. This discussion is ex- 
tended to a bilayer graphene/ (monolayer) graphene sys- 
tem in Sec. III. Corrections to the asymptotic behavior 
in the limit of large interlayer distances are considered in 
Sec. IV. A short summary concludes the manuscript. 



II. DRAG RESISTIVITY IN 
2DEG /MONOLAYER GRAPHENE SYSTEMS 

A. Model 

In this Section, we investigate a double-layer structure 
consisting of a 2DEG, located within a quantum well of 
width w, and one n-doped layer of graphene. Both elec- 
tronic systems are separated by a spacer of width d and 
embedded in a larger structure. The relative dielectric 
constants in the different regions of the structure are de- 
noted by Ki, k 2D , «2, and K3 (see Fig. 1). 

Throughout this manuscript, we consider the case 
where both layers are within the Boltzmann regime (that 
is, the regime in which the Fermi wavevector is much 



larger than the inverse mean free path) and have the 
same temperature T. Furthermore, these temperatures 
are assumed to be low, that is, 

k B T « 4 D/S , (1) 

where k B and e| D/g denote the Boltzmann constant and 
the Fermi energies of the 2DEG and graphene layers, 
respectively. In what follows, graphene is assumed to be 
the active layer, while the 2DEG is taken as the passive 
one. 

A peculiar property of the considered double-layer 
structure is the presence of interactions between two 
kinds of carriers, massive and massless fermions. Indeed, 
the carriers in the 2DEG are massive fermions with ef- 
fective mass m* and dispersion relation 

n 2 k 2 

(2) 



2m* 1 

while the carriers in the graphene layer are massless 
fermions with Fermi velocity t>§ ss 10 6 m/s and the linear 
dispersion 
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(3) 



For most practical situations, the interlayer distance 
(d) is such that the interlayer Coulomb interaction is 
weak. Thus, a lowest-order perturbation theory in the 
interlayer potential suffices and the transresistivity p 1 ^ is 
found to be given by, 5 ' 6 
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where S is the cross section area of the layers, U 2De (q, to) 
is the screened interlayer potential between the 2DEG 
and graphene layers, and <7 2D/g and r* D/g (q, ui) denote 
the Drude conductivity and the i-th component of the 
nonlinear susceptibility in the 2DEG and graphene lay- 
ers, respectively. The Drude conductivities are given by 
02D/ g = e 2 eF D/B r 2D/g /(7rfi, 2 ), where the momentum relax- 
ation times (at the Fermi energy), t 2D/s , are defined be- 
low. Furthermore, we assume that there is no electron 
tunneling between both layers, so the Fermi energies e| D/g 
can be set independently from each other in each layer. 
Due to the factor 1/sinh 2 (fiw/2fc B T), only small values 
of uj contribute to the transresistivity at low tempera- 
tures, while the screened interlayer potential (see below) 
restricts the momentum integration to small values of q. 
Therefore, we can approximate the nonlinear susceptibil- 
ities by their respective expressions in the limit of low 
energies and long wavelengths. 



B. Nonlinear susceptibilities 

Before we continue with those expressions, we briefly 
mention the general expression for the nonlinear suscep- 
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tibility within the Boltzman limit, that is, the regime of 
k F l > 1 or lot S> 1 with t, I = v f t, k F , and v F being the 
scattering time at the Fermi level, the mean free path, 
and the Fermi wavevector and velocity, respectively (for 
brevity, we suppress the index label denoting the system 
here and in the following). In this limit, the nonlinear 
susceptibilities of monolayer and bilaycr graphene as well 
as 2DEGs can be written as 



r(q,w) = 
Sh 
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(5) 



where e = |e| denotes the absolute value of the electron 
charge, A and A' band labels, eA.k = Aet the energy in 
a given band with being the dispersion of the sys- 
tem investigated [that is, Eq. (2) for 2DEGs and bilayer 
graphene and Eq. (3) for monolayer graphene] , and /* the 
Fermi-Dirac distribution function for the energy e^.k- For 
monolayer and bilaycr graphene, A = ±1 describes the 
valence and conduction bands, while for 2DEGs A = 1. 
The spin degeneracy is described by the factor g s = 2 
for 2DEGs as well as monolayer and bilayer graphene, 
whereas g v describes the valley degeneracy factor, which 
is g v = 1 in 2DEGs and g v = 2 in monolayer and bilayer 
graphene. The factor i^'k+q' which arises due to the 
overlap of the wavefunctions, is unity for 2DEGs, but 
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(7) 



for monolayer and bilayer graphene, respectively. 16 Here, 
k is the azimuthal angle of k in momentum space. We 
have also introduced the quantity /^ k k+q , which reads as 

M k ,'kVq = - [r(k)k - AAV(k + q) (k + q)] (8) 

m 

for 2DEGs and bilayer graphene and 



hv F 



r(k). 



AA'r(k + q). 



k + q| 



(9) 



for monolayer graphene. Equations (8) and (9) con- 
tain the scattering time r(k), which in general can be 
momentum-dependent. Since we are interested in the 
limit of low temperatures, we will use the expressions ob- 
tained from Eq. (5) for T — > and low energies and long 
wavelengths for the rest of the manuscript (from here on, 
we also restore the index label denoting the system). 

In 2DEGs, the main effects of both short range and 
screened Coulomb impurities can be properly described 
by considering the relaxation time t 2D to be momentum 
independent. Here, r 2D denotes the relaxation time at 



the Fermi level — a consequence of the fact that at low 
temperatures the nonlinear susceptibility is determined 
by electrons at the Fermi surface. In this case, the non- 
linear susceptibility in the limit of low energies and long 
wavelengths reads as 5,6,20 



where 



T 2D (q,w) = 
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and v F D and k F D are the Fermi velocity and wavevector 
in the 2DEG layer. 

Contrary to the case of the 2DEG, the relaxation 
time describing electron-impurity scattering in graphene, 
which is proportional to the momentum, that is, r g (k) = 
Tofc, with To being a constant of proportionality, 15,21-25 
is widely used as a model for the relaxation time in 
graphene. In this case, the nonlinear susceptibility in 
the limit of low energies can be written as 17-19 



r* (q,w) 



4euJTg q l 
hmv% q 
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where 



i- 



6 (2fc| - q) 



(13) 



Here, v§, and fcf denote the Fermi velocity and wavevec- 
tor in graphene and the relaxation time at the Fermi level 
is given by r g = r fc|. We note that Eqs. (12) and (13) 
are the same results one would have obtained if the re- 
laxation time in graphene had been assumed as constant 
(that is, r g (k) = r g = const.). 12 Indeed, it has been 
shown in Refs. 17-19 that — if isotropic relaxation times 
are assumed — the form of r g (|k|) as a function of the mo- 
mentum does not affect the low-temperature limit of the 
nonlinear susceptibility, Eq. (12), and one can replace the 
momentum-dependent relaxation time by its value at the 
Fermi level. Moreover, one can notice that, within the 
limit of Eqs. (10) and (12), the momentum integration is 
cut off for q > 2kT ,s . 



C. Interlayer potential and transresistivity 

The screened interlayer potential can be found by solv- 
ing the corresponding Dyson equation and can be written 
as 



^2 Dg (q,w) 



U^i (q) 
s 2 D g (q,w) 



(14) 
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with 

e 2Dg (q, u) = 1 + U$ (q) II 2D (q, w) 

x [l + C/ g <o) (q)n s (q,w)] 

2 

c^Dg (q) n 20 (q, u) n g (q, w) . 



(15) 



where n 2D / g (q, w) are the polarization functions of the 
individual layers, for each of which we use the respec- 
tive expressions obtained from the random phase ap- 
proximation (RPA) at zero temperature. The bare in- 
tralayer and interlaycr Coulomb potentials can be writ- 



ten as Ui°' /S (q) = (47re 2 /g)/ 2D/g (qd, qw) and U^i (q) = 
(8ne 2 /q)f 2Dg (qd, qw), respectively, where the form fac- where Q 



(0) 



tors / 2D / g (qd, qw) and / 2Dg (qd, qw) arc determined by 
solving the Poisson equation of the system (see Ap- 
pendix A). Since at low temperatures only small values 



of u) contribute to Eq. (4), we approximate the dynamic 
by the static polarization functions, that is, we replace 
U 2Dg (q, w) by the static interlayer potential U 2Dg (q, 0). 

Within the above approximations, the transresistivity 
is diagonal because the nonlinear susceptibilities and the 
screened interlayer potential are isotropic, and it is there- 
fore enough to calculate p D = p* x . Thus, we obtain the 
transresistivity at low temperatures, 



h 4tt (k B T) z F (Q™d, Q% F d, w/d) 
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and 



denote the bare 



Thomas-Fermi wavevector and the (total) density of 
states at the Fermi level in each individual layer. The 
function F(y 2T> ,y s ,r) is given by the integral 
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where y 2U / s = Q 2 ° F /s d, r = w/d, x = qd, and 
n 2D (a;/d,0) 
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l + e 



i 



(- 

\2k F d J Ak F d 
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(19) 

are the static, dimensionless polarization functions of the 
2DEG 20 ' 26 and graphene, 27 ' 28 respectively. 



D. Asymptotic behavior 

In general, the integration in Eq. (17), and therefore 
the transresistivity in Eq. (16), have to be computed nu- 
merically. However, simplified analytical expressions de- 
scribing the asymptotic behavior of the drag resistivity 
as described by Eqs. (16) and (17) can be obtained in the 
limits of large and small interlayer distances by replac- 
ing each of the different relative dielectric constants K\, 



K2, K3, and k 2D by an average relative dielectric constant 
of the entire structure, k. Then, we can introduce the 
screened Thomas-Fermi wavevectors, <?tf /s = <9tf /e /k- 
Moreover, the form factors reduce to 



/2D — 
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(20) 
(21) 

(22) 



Since the upper boundary of the integral given by 
Eq. (17) is restricted by the minimum of the Fermi 
wavevectors [min (2k F D d, 2k F d)], the polarization func- 
tions can be replaced by their long wavelength limits, 
that is, n 2D/g (x) — > 1. 

Below, we study the asymptotic behavior of Eqs. (16) 
and (17) in three different limits corresponding to small 
and large values of q^ /s d. Here, we note that Eqs. (16) 
and (17) have been derived under the assumption that 
the nonlinear susceptibilities can be approximated by 
Eqs. (10) and (12). As shown in Rcf. 17, however, this is 
not the case for the limit d — 0, where the Fermi energy 
is no longer the largest scale of the system and Eq. (12) 



5 



system 


Pd 






(active-passive) 


(q&d < 1, K /P d < 1) 


(<£ F p d< 1, fc F /p d> 1) 


(<g£d> 1, K /P d> 1) 


massive-massive 


~ 1 

"I 


ln[(^ F +9^. F )d] 
(n a «„)3/2 




massless-massless 


oc — 

1 a 


ln K%A T 7+A/' 1 p) ti / K l 

OC ; 




massless-massive 


oc 4, 


oc v r ^- TF 

V a p 


oc l/[Kn f ) 3/J «i 4 ] 



TABLE I: Asymptotic behavior of the transresistivity p D as a function of the densities and interlayer distance for different 
systems. In the limit q^! F d 2> 1, the three systems exhibit identical behavior in the leading order and one needs to consider the 
subleading correction Ap D in order to see differences in the transresistivity (see Sec. IV). The subscript a (p) refers to the active 
(passive) layer. The screened Thomas-Fermi wavevectors, particle densities, interlayer distance, and average dielectric constant 
are denoted, respectively, by q^ F , n a / p i d, and k. We have also introduced the constant a — 4:e 2 /(Khv F ) with v F denoting the 
Fermi velocity of the massless particles. Since there is no general analytical formula for the small interlayer separation limit, 
q^ F d <C 1 and k F /v d <C 1 (first column), we provide formulae for the case of high densities and k F = fep, that is, n a = 2n p for 
a 2DEG-graphene system and n a = n p for a BLG-graphene system (see Appendix B). 



is not a good approximation for the nonlinear suscepti- 
bility in graphcne [the same is also true for the nonlinear 
susceptibility of the 2DEG given by Eq. (10)]. Only for 
weak interaction strength, Eq. (12) can describe the non- 
linear susceptibility for small, but finite d reasonably well 
[For d = 0, the transformation in the integral of Eq. (17) 
cannot be used.]. 18,19 Thus, the limit of small interlayer 
separation presented in the following should be under- 
stood in this way. 



1. 



Small interlayer separation limit (q^, F /g d, fc F D d <C 1) 



In this case, the integration in Eq. (17) is restricted 
by the upper boundary xq = min (2k F °d, 2k F d) and we 
obtain 



p D = 



h (k B T) q™q% F n 



a 1 £ 2D^S jU2D L>& lO 
C C F C F "-p n>F 



/ (fc F D , q™ + <£ p ) - g (fcf, fc FJ q™, c£ P ) ^kf^w + O (k?k*w 2 ) , (23) 
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(24) 



(25) 



with j/o = min [y/kf/k%, y/k%/kfj and r 2D/g = 

<7tf /b / \/k' F D k F . Equation (23) shows that in the small 
interlayer separation limit the transresistivity does not 
depend on d. Such a behavior has also been found 
in graphene-graphene double-layer structures. 18,19 More- 
over, g (fc F D , fc|, g™! 9tf) > an d thus the transresistivity 
is reduced for finite widths of the quantum well. Equa- 
tions (24) and (25) have — in general — to be computed 
numerically. However, for certain limiting cases analyt- 
ical formulas can be derived for which we refer to Ap- 



pendix B. In particular, when the particle densities are 
such that n = n s = 2n 2D we obtain, in the leading or- 
der of 1/n, p D cx 1/n 2 . This dependence can be seen as 
an intermediate behavior when compared to the results 
expected for 2DEG-2DEG (p D oc 1/n 3 ) and graphene- 
graphene (p D cx 1/n) double-layer structures when the 
particle density is equal in both layers. 17-19 
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2. Intermediate limit (q™ /s d <1< fcp D/g d) 

Whereas the limit considered above requires small 
interlayer distances, d must not be too small for the 
limit q^ /s d <C l,kl D/s d 3> 1. It is difficult, however, 
to reach this limit experimentally because in graphcne 
oc fcf /k. Thus, this limit can only be reached if k is 
very large. 

Since kp D/e d ^> 1, the integral in Eq. (17) is practically 
restricted by the Coulomb interaction with the main con- 
tribution arising from values x = qd < 1. Therefore, we 
can approximate A 2D/g (a;/d) by A 2D/ls (x/d) — > 1 and find 



h (k B T) z q™q* 



hi 



fcjPfcf 12 



1 



-J + o [ Wrf) 2 ] 



(<g° + <?f F )d 

(26) 

which has been expanded in powers of w/d. From 
Eq. (26), we obtain p a ~ ln[(ay / 7rrl^ + q™)d]/y/ri$ 
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for the dependence of p D on the carrier densities, with 
a = Ae 2 /{nhv%). 
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FIG. 2: (color online). Dependence of the transresistivity on 
the temperature, T, for air/graphene/A^Oa/GaAs/AlGaAs 
and air/graphene/Si02/GaAs/AlGaAs structures with an in- 
terlayer distance of d = 20 nm for different widths of the GaAs 
quantum well (w = 0,5,10 nm). The electronic densities of 
graphene and GaAs have been set to the values n g = 1.0 x 10 12 
cm -2 and riGaAs = 1.0 x 10 11 cm -2 , respectively. 



E. Numerical calculations 



3. Large interlayer separation limit (g^p /E d, fcp D/s d 3> 1) 

As for the limit of intermediate interlayer separation 
above, the main contribution to Eq. (17) arises for x < 1 
and we can approximate A 2U / s (x/d) — ¥ 1. In this case, 
the values contributing to Eq. (17) satisfy x <C gto/g/ft 
and we obtain 

h_ (k B T) 2 ttC(3) 

Pd e 2 e™e s F (k™d){k e F d)(q™d)(q s TF d) 32 

J (720tt 2 + 1350) C(3) - 7T 4 (4tt 2 - 15) w 
X j 1 540tt 2 C(3) d 

+ 0[(w/df]^. 

(27) 

Consequently, the dependence of p D on the carrier densi- 
ties is given by p D ~ l/[(ri 2D n g ) 3 / 2 d 4 ], which is the same 
asymptotic behavior as one would expect in a double- 
layer structure consisting of two 2DEGs or one consisting 
of two graphene layers. Only when higher order terms in 
the series expansion of A 2D/g (a;/d) [see Eqs. (11) and (13)] 
are taken into account, one can find a difference in the 
asymptotic behavior (see Sec. IV). 

The asymptotic behavior of the transresistivity as a 
function of the carrier densities is summarized in Table 
I, where, for comparison, the results corresponding to 
massive-massive and masslcss-masslcss systems have also 
been included. 



We have performed numerical calculations us- 
ing Eqs. (16) and (17) for two different struc- 
tures, air/graphcne/A^Oa/GaAs/ AlGaAs and 
air/graphene/Si02/GaAs/ AlGaAs, in which graphene 
plays the role of the active layer and the 2DEG formed 
in the GaAs quantum well constitutes the passive one. 
The two structures differ in the materials conforming the 
spacer, AI2O3 in the former and Si02 in the later case, 
which possess different dielectric constants, K2 = 9.1 
and k 2 = 3.9, respectively. The remaining system 
parameters used in the evaluation of the transresistivity 
are K\ = 12.9, k 3 = 1, « 2D = 12.9, and the electron 
effective mass in GaAs, rn* Cahs = 0.063 mo. Here, mo 
represents the bare electron mass. 

The temperature dependence of the transresistivity in 
the two considered structures is shown in Fig. 2 for dif- 
ferent widths of the GaAs quantum well. The interlayer 
distance is d = 20 nm and the densities nc a As = l.Ox 10 11 
enr 2 and n„ = 1.0 x 10 12 cm -2 . Here, the densi- 
ties chosen for each layer reflect the possibility that for 
double-layer systems consisting of 2DEGs and graphene 
one might be able to have a combination of rather differ- 
ent densities in both layers, a case difficult to achieve if 
both layers consist of the same material. As can be seen 
from Eq. (16) and Fig. 2, — p D ~ T 2 . This is the same 
temperature-dependence found for the transresistivity in 
ballistic 2DEG-2DEG bilayers. 2-6 One can also appreci- 
ate in Fig. 2 that, for a given temperature, the smaller the 
well width, the larger the size of the resistivity. This is a 
general behavior, which, according to Eqs. (23), (26), and 
(27), occurs in both the q^ /e d < 1 and the q™' & d > 1 
limits. Note, however, that for the set of parameters 
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FIG. 3: (color online). Dependence of the transre- 
sistivity on the electronic density in graphene, n s , at 
T — 100 K for air/graphene/Al 2 3 /GaAs/AlGaAs and 
air/graphene/Si02 /GaAs/ AlGaAs structures with an inter- 
layer distance of d — 20 nm and for different widths of the 
GaAs quantum well (w = 0, 5, 10 nm). The electronic density 
of GaAs has been set to the value nGaAs = 1.0 x 10 11 cm -2 . 



III. DRAG RESISTIVITY IN BILAYER 
GRAPHENE/MONOLAYER GRAPHENE 
SYSTEMS 



Apart from the 2DEG-graphene system considered in 
the previous section, the Coulomb drag between mass- 
less and massive fermions may also be realized in a 
double-layer structure consisting of bilayer and mono- 
layer graphene isolated from each other by a spacer. 
Compared to the system investigated in Sec. II, the quan- 
tum well of width w containing the 2DEG is replaced by 
a sheet (w = 0) of bilayer graphene. As before, mono- 
layer graphene is assumed to be the active layer, while 
bilayer graphene is taken as the passive one. Likewise, 
the ballistic case is assumed for both layers and only low 
temperatures are considered. We restrict our analysis to 
the case in which both layers are electron doped. In such 
a case, bilayer graphene consists of an electron band with 
a parabolic dispersion as in Eq. (2), but with the BLG 
effective mass m£ g = 0.033 mo. With the relaxation time 
in bilayer graphene at the Fermi level given by T bg , for 
T — > one obtains 



r « «■">- -St ?**<•> (28) 



considered in Fig. 2, which corresponds more to the limit 
q^ /e d 3> 1, both — p D and its changes with the well width 
are larger in the air/graphenc/AkjOa/GaAs/ AlGaAs sys- 
tem than in air/graphenc/Si02/GaAs/AlGaAs. This can 
be qualitatively understood by noting that in such a 
limit both —po and —dpo/dw are proportional to 7c 2 
[see Eq. (27), where ~k enters via the screened Thomas- 
Fermi wavevectors] . Consequently, since n-i (and there- 
fore k) is larger in air/graphene/Al 2 03/GaAs/AlGaAs 
than in air/graphene/Si02/GaAs/AlGaAs, both the ab- 
solute size of the transresistivity and its changes with w 
are expected to be larger in the former structure com- 
pared to the later one, as is indeed seen in Fig. 2. 



The dependence of the transresistivity on the den- 
sity in the graphene layer at T = 100 K is 
shown in Fig. 3 for graphene/ AI2O3 /GaAs/ AlGaAs and 
graphene/Si02/GaAs/AlGaAs structures. Here, the in- 
tcrlaycr distance is set at d = 20 nm and the density in 
the GaAs layer at riQaAs = 1-0 x 10 11 cm -2 . We can 
fit the curves in Fig. 3 to —po oc ng and extract values 
between j3 —1.19 and j3 f=s —1.30, which are closer to 
—3/2 than to —1/2. This is consistent with the fact that, 
for the parameters chosen, we are approaching the limit 
of fe /g rf> I- 

Regarding the dependence of — p D on the well width, 
both Fig. 2 and Fig. 3 exhibit the same qualitative behav- 
ior: the smaller the well width, the larger the absolute 
value of the transresistivity. 



for the nonlinear susceptibility in the limit of low energies 
and long wavelengths, where 



A bg (g) 



X 6 (2k** - q) 




(29) 

and fcp E and i>p E denote the Fermi wavevector and velocity 
in bilayer graphene. 

We can replace the quantities describing the 2DEG 
layer by the respective quantities of bilayer graphene and 
use the results from Sec. II. This means that, aside from 
setting w = 0, replacing the Fermi energy, the Fermi 
velocity, the Fermi and Thomas-Fcrmi wavevectors, as 
well as A 2D , the polarization function of the 2DEG in 
Eq. (17) has to be replaced by the polarization function 
of bilayer graphene calculated in Ref. 29. However, in 
contrast to a 2DEG, one has to take into account the 
valley degeneracy of bilayer graphene, the net effect of 
which is an additional factor of 1/2 in Eqs. (16), (26), 
and (27). 

Whereas the limiting cases of intermediate and large 
interlayer separation from Sec. II D do not depend on 
the exact form of A 2D/bg because of kl D/b9 d » 1 and can 
therefore be described by Eqs. (26) and (27), respectively 
(and taking into account the additional factor of 1/2 for 
bilayer graphene-graphene systems as well as setting uu = 
0) , this cannot be done in the small interlayer separation 
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FIG. 4: (color online), (a) Dependence of the transresistiv- 
ity of an air/graphene/AbOa /bilayer graphene/Si02 struc- 
ture on the electronic density of graphene (n g ), for different 
electronic densities in bilayer graphene (nb g ) at T = 100 K. 
(b) Temperature dependence of the transresistivity for the 
same structure as in (a) at fixed graphene electronic density, 
n E = 5x 10 11 cm~ 2 . In both cases, (a) and (b), the interlayer 
distance is d = 20 nm. 



limit. For this limit, we obtain 



h (k B T) 2 q*%q* F ir 



e 2 e F s e F k F e ki 24 



f {k F , k F , q tf ~\~ ^xf) 



(30) 



where f(k F e , k F , g TF ) is given by Eq. (24) and 
fi (k F s j k F , 9tf) = 

4fc| r yVg - (fc^/fc F )y 2 ] [1 - (g/fcgjgj 
r y + g TF /(2^ i fcf) 



(31) 



with yo = mm ( V fc F E /fc|, VWW) . In general, Eq. (31) 
has to be computed numerically, but it is possible to de- 
rive analytical formulas for certain limiting cases (see Ap- 
pendix B). 

Thus, the transresistivity due to the Coulomb drag be- 
tween massive and massless fermions in 2DEG-graphene 
and BLG-MLG structures is characterized by the same 
generic expressions for the intermediate and large inter- 
layer separation limits. For small interlayer distances on 
the other hand, there is an additional contribution in 
the case of the BLG-graphene system compared to the 
2DEG-graphcne system [see the massless-massive system 
in Tab. I]. 

Figure 4 shows the dependence of the transresistivity 
in an air/graphene/A^Oa/bilayer graphene/Si02 struc- 
ture on the electronic density of graphene (n g ) at T = 100 
K [Fig. 4 (a)] as well as on the temperature [Fig. 4 (b)] 
for n s — 5 x 10 11 cm -2 and different densities in bilayer 
graphene (n bs ). As anticipated above, the qualitative 



trends displayed in Fig. 4 are similar to those shown in 
Figs. 2 and 3. However, for the set parameters taken in 
Fig. 4, the system starts to approach the large interlayer 
separation limit and the absolute value of the drag resis- 
tivity in the air/graphene/A^Os/bilayer graphene/Si02 
appears to be smaller than in the 2DEG-graphcne sys- 
tem. This behavior can be understood from Eq. (27) 
by taking into account that the screened Thomas-Fermi 
wavevector in BLG is larger than in GaAs. 



IV. CORRECTIONS TO THE LARGE 
INTERLAYER SEPARATION LIMIT 

As mentioned in Sec. II and shown in Tab. I, the 
asymptotic behavior of the transresistivity in a massless- 
massive double-layer system in the limit of q^J F d 3> 
1 is identical to the behavior in massive-massive and 
massless-massless systems [a/p denote the active/pasivc 
layers). Only when higher order terms in the series 
expansion of A a/p (x/d) [see Eqs. (11), (13), and (29)] 
are taken into account, one can find a difference in the 
asymptotic behavior. 

In general, we find that the leading correction (for w = 
0) to Eq. (27) is given by 



h 5ttC(5) (k B Tf [c p (k F d) 2 + c.(fcgrf) 2 ] 
e 2 256g a g p e F e F (fc F ci) 3 (fc F d) 3 (q^ F d) (q? F <f) 



(32) 

where a/p, fc F /p , and q F /p denote the active/passive layers 
and their respective Fermi and screened Thomas-Fermi 
wavevectors. The parameters c a/p and g a/p are specific 
of the system comprising the active/passive layers, with 
c 2D = 1, Cg = —1, and c bg = —7, as well as g 2 u = g s = 1 
and <7b g = 2. 

From Eq. (32), we find Ap D ~ (n s — 
2"-2D)/[("-2D"- g ) 5 ^ 2 ^ 6 ] for the density depen- 
dence in 2DEG-graphcnc systems in contrast to 
Apo - =F(n* + n p )/[(n a n p ) 5 / 2 >] in 2DEG-2DEG (-) 
and graphene-graphenc (+) systems. In particular, 
under the condition of n = n & = 2n 2D the correction 
to the drag vanishes for the 2DEG-graphene system 
but remains finite, with the asymptotic behavior 
A Po - l/(n 4 d 6 ), for the 2DEG-2DEG and graphene- 
graphene systems. Similarly, Eq. (32) can be used to 
describe deviations from the asymptotic behavior for 
2DEG-BLG, BLG-BLG, and BLG- (monolayer) graphene 
systems. 

In the limit q^pd 3> 1, the drag correction is small, in 
agreement with the trend described by Eq. (32). For the 
set of parameters considered here (which correspond to 
a region close, but still not in such a limit), we have 
found from our numerical calculations that while the 
drag correction in air/graphene/Al 2 03/GaAs/AlGaAs 
and air/graphene/Si02/GaAs/AlGaAs turns out to be 
still small (a few percent of the total drag), it becomes 
relevant for the case of air/graphene/A^Oa/bilayer 
graphcnc/Si02, in which it represents about 30% of the 
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total transrcsistivity. 

For the case of air/graphene/Al 2 03/bilayer 
graphenc/Si02 with the parameters used in Fig. 4, 
the limit q^^d 3> 1 is reached when the interlayer dis- 
tance is increased to values d > 80 nm. In such a limit, 
we found a very good agreement between our numerical 
calculations and Eq. (32). Our calculations indicate that 
the drag correction decreases from 12 — 13% to 1% (at 
n g = 1.5 x 10 12 cm -2 ) of the total transrcsistivity when 
the interlayer distance is increased from d = 20 nm to 
d = 80 nm. 



V. CONCLUSIONS 

In this manuscript, we have studied the Coulomb drag 
at low temperatures in a double-layer structure composed 
of a 2DEG and a graphene layer, both of which were 
treated as being in the Boltzmann regime. We have 
written down a formula to describe the transresistivity 
of such a system at low temperatures and have analyzed 
the temperature and density dependence of this formula 
both analytically and numerically. Analytical formulae 
have been derived to describe the asymptotic behavior 
in both, the small and large interlayer separation lim- 
its and compared to the respective behavior in massive- 
massive as well as massless-massless systems. It has been 
found that for q^^d <C 1 each system, massive-massive, 
massless-massless, and massless-massive, possesses a dif- 
ferent dependence on the carrier densities, whereas the 
three systems share the same behavior in the dominant 
contribution to p D for q^d 3> 1. Only looking at higher- 
order corrections allows us to distinguish between the 
different systems in this regime. Furthermore, the ef- 
fect of a finite width of the quantum well in which the 
2DEG is formed has been investigated and we have seen 
that with increasing well width the absolute value of the 
transresistivity is reduced. Finally, we have also studied 
a BLG-graphcne system, which we found to be qualita- 
tively similar to a 2DEG-graphene system in the large 
interlayer separation limit, but different in the limit of 
small interlayer separation. 
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Appendix A: Bare Coulomb potential 

The bare Coulomb potentials can be obtained from the 
Poisson equation, which in cylindrical coordinates (p, z) 
reads as 

V [«(z)V0(p - p'; z, z')] = AveS (p - p') 5 (z - z) 

(Al) 



for a point charge located in a geometry as shown in 
Fig. 1. Here, the relative dielectric constant is given by 



k(z) 



/*3 for z > d + w 

K2 for w < z < d + w 

k 2D for < z < w 

K\ for z < 0. 



(A2) 



Introducing the Fourier transform of <p with respect to the 
in-plane coordinates p and insertion in Eq. (Al) yields 



d 



k{z) 



d0(q;z,z') 



dz 



■ n(z)q 2 (f>(q; z, z') — 4ire5 (z — z') . 

(A3) 

This equation is solved for each region given in Eq. (A2) 
(and each combination of z and z') separately and we re- 
quire the global solution to be continuous and its deriva- 
tive to be piecewise continuous with a jump of Aire at 
z = z' . 

Having determined the potential 0(q; z, z') in this way, 
the bare Coulomb potential can be calculated from 

00 00 

C/f (q) = -e J dz J dz'^z,z')\ Xt {z)\ 2 \ X3 {z')\ 2 , 

— OO —OO 

. (A4) 

where Xi/j( z ) describes the localization in z direction of 
a particle located in the 2DEG (i = 2D) or graphene 
(i = g) layers. For graphene, we assume the electrons to 
be perfectly localized and therefore 



IXgWI 2 = S(z-d-w). 



(A5) 



The transversal wavefunction of an electron located in 
the 2DEG quantum well, on the other hand, is assumed 
to be given by that of the ground state of an infinite 
one-dimensional potential well, 



|x»d(«)I = - shl {— 6(w-z)Q(z), (A6) 
w \ w J 

that is, we assume that only the lowest quantum well 
subband is occupied. From the solution of Eq. (A3), 
0(q; z, z'), and Eqs. (A4)-(A6) we find the interlayer po- 
tential f/ 2 Dg (q) = (87re 2 /g)/ 2Dg (qd,qw) (i = 2D and j = 
g or vice versa) and the intralayer potentials J7 2 D/g (q) 



(47re 2 /g)/ 2D/g (qd,qw) (i = j 
the form factors 



2D and i = j = g) with 



/2D g (x, y) 



/ g {x, y) = 



27T 2 K 2 {«i [cosh(y) - 1] + k 2D sinh(y)} 



y(y 2 +4n 2 )Af(x,y) 



K2 cosh(x) [ki sinh(y) + k 2D cosh(y)] 



(A7) 



+ 



Af(x, y) 

k 2D sinh(a;) [ki cosh(y) + k 2D sinh(y)] 



Af(x,y) 



(A8) 
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/sd {x, y) 



where 



K X K 2 [«2 sinh(a;) + k 3 cosh(a;)] {64tt 4 [1 - cosh(y)] + y (y 2 + 4tt 2 ) (3y 2 + 8?r 2 ) sinh(y)} 

2n 2D y 2 (y 2 +^) 2 M(x,y) 
[k 2 («i + k 3 ) cosh(x) + («| + K!« 3 ) sinh(a;)] [y (32tt 4 + 20ir 2 y 2 + 3y 4 ) cosh(y) - 32tt 4 sinh(y)] 

2 y 2 ( y 2 + 4 7r 2) 2 ^ a . )J/ ) 

k 2D [k 2 cosh(a;) + k 3 sinh(x)] y (y 2 + 47r 2 ) (3y 2 + 87r 2 ) sinh(y) 
+ 2y 2 (y 2 + 4ir 2 ) 2 Af(x,y) ' 



Af(x, y) =k 2 cosh(x) [k 2D (k x + k 3 ) cosh(y) + («i« 3 + k 2D ) sinh(y)] 

+ sinh(x) [k 2D (k 2 + «i/t 3 ) cosh(y) + (kiK 2 + k 3 k 2 d ) sinh(y)] 

I 



(A9) 



(A10) 



In the limit of y — >• [/i(x) = /i(x, 0)], which describes 
the setup of Sec. Ill, we recover 



K2C 



(«i - k 2 ) (k 2 - 


- «s) 


f e 22 («i 4 


" K2) («2 


+ «a)' 
(All) 


(«2 - 




- («! + K 2 ) 


e 2 * 




(«1 - «2) («2 - 


k 3 )4 


■ e 2x (ki + 


K2) («2 ^ 


-«3)' 

(A12) 


(«2 " 


- «3) " 


\~ («2 + «3 


e 2x 




(«1 - K 2 ) («2 - 


- «3) " 


h e 2x (ki + 


K 2 ) («2 " 


f «a)' 
(A13) 



/2D g (a;) 



/ 2D (x) 



consistent with the results found in Refs. 30,31. 



Appendix B: Limiting cases for small interlayer 
distances 



As noted in Sees. II and III, the integrals defined by 
Eqs. (24), (25), and (31) have to be computed numeri- 
cally, but it is possible to obtain analytical formulas for 
certain limiting cases. In the following, some of these 
cases are presented. 

If k F w k F = k F (m = 2D/bg), that is, for carrier 
densities n s 2n 2D in a 2DEG-graphene structure or 
n s « n bg in a BLG-graphcne structure, the integrands 
can be expanded around k F /k F = 1 and the integration 
of the lowest order contributions yields 



/ (fc F , k F ,q TF ) = In 



1 + 7 



1 



and 



fi{k F ,k F ,q TT ) =4 7 2 (3 - 5 7 2 )ln 



1+7 
l + 7 " 



(Bl) 



7 
16 



(B2) 



+ 207 J - IO7 2 - — 7 + 1 



where 7 = q TF /(2k F ). For the massive-massless systems 
investigated in Sees. II and III, the functions / and /1 



have to be evaluated at q TF = g« p + g™ F [see Eqs. (23) 
and (30)] and thus 7 = a/2 + q^ F /(2 y /nn^). Assuming 
that we are in the high-density limit (consistent with our 
usage of RPA), we can further simplify the transresistiv- 
ity, given by Eqs. (23), (30), (Bl) and (B2) if kg = fcf, 
by expanding it in powers of the density, from which we 
obtain p D oc 1/n 2 . This is the formula used in Tab. I for 
massive-massless systems. 

In the case of kg /k% <C 1, we can expand the inte- 
grands around kg /k F = and then perform the integra- 
tion of the lowest order contributions analytically, 



/ (k F , k F , q TF ) 



1 e ( i-7<)m(i±^S) 



e ( 7 < - 1) 



| 7 |-1| 3/2 



arctan 



|7<-1| 3/2 



(B3) 



and 



fi (K,k F ,q TF ) = 
2 7 2 (47 2 -3) 



+ 29( 7 <-l) 



6(l- 7< )ln 



1- vgEg 

1+ n/T+tT 



arctan 



1 



4 



(B4) 



+ 8tt 7 < - I67 2 - 2tt7< + 



where 7< = q TF /(2kg). 

Finally, we consider the case of kg/k F ^> 1. Here, we 
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can expand the integrands around /k™ = and obtain 



/ (fc F j fcp, q TF ) — 



where 7> = q TF /(2k%). Moreover, the main contribu- 
tion in this limit arises from / because / (fc™, fcg, q TF ) ^> 



2-) 



-9(l- 7 >)ln 



i + Vi 



7^ 



7> 



+ e ( 7 > - 1) 



arctan 



(B5) 



+ ^7> ~ 2, 
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